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We study dimuonium {fi'^ production in the quark-gluon plasma created in relativistic heavy 
ion collisions. The production is controlled by the process qq — ^ {ii^n~)g, and the dimuonium 
motion in the plasma is described by a transport equation. While the electrodynamics dominated 
dimuonium yield is not high enough, the transverse energy distribution carries the information of 
the plasma at RHIC and LHC energies. 
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PACS numbers: 36.10.Ee, 25.75.Cj, 12.38.Mh 



The observation of positronium (e+e ) 1] and the clar- 



ification of the leptonic nature of muons lead to a natural 
question that whether there are "true muonium" [2] or 
"dimuonium" [3, 4] states (fi^n^). Note that the termi- 
nology "muonium" has been signed to the bound state 
(/i+e"), which has been discovered in 1960 Since 
dimuonium is an ideal system of quantum electrody- 
namics (QED), many of its properties are theoretically 
predicted [1, |g| with standard calculations of QED. For 
instance, dimuonium is one of the most compact QED 
systems with Bohr radius 512 fm, the binding energy 
of the ground state (Is) is 1.4 KeV, and the life times 
for the states ^5*0 and ^Si are respectively 0.602 ps (de- 
cay to 77) and 1.81 ps (decay to e+e"), which are both 
much shorter than the life time of muon (2.197 /is) and 
thus the muon weak decay can be ignored in the pro- 
duction of dimuonium. Many dimuonium production 
mechanisms have been proposed, like direct fJ-'^ col- 
hsion I2j], e+e" ^ (a^+M") [7], t^'P ^ H, 
-yA {^+^^)A [6] and eA -J> e(/i+^")A g, where 
A stands for a nucleus. Dimuoniums can also be pro- 
duced in heavy ion collisions through a pure electromag- 
netic process A1A2 — AiA2(/x+/i~), where the mecha- 
nism is the fusion of coherent photons emitted from the 
nuclei Q. Recently, the possibility of dimuonium pro- 
duction at modern electron-positron colliders is investi- 
gated [loj . However, dimuonium has not yet been exper- 
imentally discovered. 

It is generally believed that there is a quantum chro- 
modynamics (QCD) phase transition in hot and dense 
nuclear matter, which is related to the deconfinement 
process in movirig from a hadron gas to a quark-gluon 
plasma (QGP) |11[. The realization of such a phase 
transition in laboratories can only be through relativis- 
tic heavy ion collisions [12]. Considering the decreasing 
temperature and density during the rapid expansion of 
the fireball formed in the collisions, the QGP can not 
be measured directly in the cold and dilute final state, 
and one needs signatures to identify the QGP formation. 
Since leptons interact with particles only through elec- 
tromagnetic channel, the leptons produced in the QGP 
carry the information of the plasma and will not lose it 
when they pass through the collision region to reach the 
detectors. Taking into account the fact that the lepton 
number in the fireball increases with the temperature of 



the plasma, the probability to produce a dimuonium state 
in high energy nuclear collisions should be much larger 
than that in elementary electron-positron and nucleon- 
nucleon collisions. In this paper we study the dimuonium 
production in a hot QGP and calculate the production 
rate in high energy nuclear collisions at relativistic h eavy 
ion collider (RHIC) and large hadron collider (LHC) 
Since the dimuonium life time is much longer than the 
QGP life time (about several fm/c), the dimuonium de- 
cay and dissociation in the QGP can safely be neglected. 
In our calculation we will also ignore the dimuonium pro- 
duction in the initial state and in the hadron gas, since 
their contribution should be much smaller in comparison 
with the production in the QGP. 

There are two kinds of production processes in QGP 
for the dimuonium states (/i^/i~), one is from the lep- 
ton pair annihilation and the other is from the quark- 
antiquark pair annihilation. While the QCD processes 
in vacuum at RHIC and LHC energies may still be non- 
perturbative, the medium effects will largely reduce the 
effective coupling constant as (^ 0.3) jl^ and may 
make the perturbative analysis available. By counting 
the number of vertexes and considering the phase space 
suppression, one can simply estimate the order of the 
dimuonium production cross sections in different pro- 
cesses. The cross sections of the two processes at leading 
order, l^l^ — > (fi'^fi") and qq — )■ {p'^fi~), are of the 
order of a^, with a being the electromagnetic coupling 
constant. Due to the requirement of energy conservation, 
only the initial states in a narrow energy window around 
the binding energy 1.4 KeV contribute to the production 
processes, and therefore the both cross sections at lead- 
ing order become very small. The next leading process is 
qq — >■ {fi~^ n~ ) g with the cross section of the order of ctgCt^ ■ 
The subsequent processes are l^l^(qq) (/i+/i^)7 and 
qq — > {p.^ ^~)gg. The former is of the order of and the 
latter is of the order of a^a^. Since the latter includes 
three particles in the final state, the process is strongly 
suppressed in comparison with the former. Other pro- 
cesses are higher order contributions. Therefore, we con- 
sider in the following the main process qq — {fi^^~)g for 
the dimuonium production in the QGP. The contribution 
from tlic other processes is at least ten times smaller from 
the simple vertex counting. 

The Feynman diagrams at tree level for the main pro- 
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duction process qq — >■ {fj,'^ fi~)g are shown in Fig[TJ Since 
strange quarks are much heavier than hght quarks, the 
amount of strange quarks in QGP is less than 10% of the 
amount of hght quarks at RHIC and LHC energies, and 
we consider only light quarks in the calculation (Nf — 2), 
and their current mass is set to be zero. The scattering 
amplitude A4{qiq2 — > PiP2k) for the process qq 
can be written as 
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where qi,q2,pi, P2 and k are respectively the momenta 
of q, q, /i"*", /i~ and g, Q is the quark charge number, t 
is the color matrix, g is related to the effective coupling 
constant via the definition as — g'^/{ATr), u,v,u and v 
are Dirac spinors, and is the gluon polarization vector 
satisfying e*(fc)e^(fc) = —g^i, + k^k^/rn?g with thermal 
gluon mass vig — 2gT /?> in the deconfinement phase 
For simplicity, we have dropped in ([T]) the flavor, color 
and spin indices of quarks and gluons. 





FIG. 1: The Feynman diagrams at tree level for the main 
dimuonium production process qq — >■ (fi'^ ii~)g in QGP. 

Introducing the total and relative momenta p = pi +p2 
and p' = (pi — p2)/2 of /Lt+ and the amplitude 
■M{qiq2 PiP2k) can be expressed as J^{qiq2 — > pp'k). 
Since the bound state {jX^ is non-relativistic, we can 
compute its production amplitude M.((iiq2 — >■ pk) by in- 
tegrating out the relative momentum p' in the center-of- 
mass frame of the dimuonium llSl , 
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where rud is the dimuonium mass, and -(ACp') is the rel- 
ative wave function for the ground state in momentum 
space with ao being the Bohr radius. From the known 
scattering amplitude, we can calculate the dimuonium 

production cross section cr^^ ^ {s) as a function of 



W. 
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52) and the transition probability 



Ap-'^i^ )g 



(s). 

The medium created in high-energy nuclear collisions 
evolves dynamically. In order to extract information 
about the medium by analyzing the dimuonium distri- 
butions, both the hot and dense medium and the dimuo- 
nium production process must be treated dynamically. 
In this paper, we treat continuous dimuonium production 
in QGP self-consistently, including hydrodynamic evolu- 
tion of the QGP. Since dimuonium is a pure electromag- 
netic system, it can not be thermalized with the medium 
which is governed by strong interactions. Thus its phase 
space distribution should be controlled by a transport 
equation. The transport equation should then be solved 
together with the hydrodynamic equation which charac- 
terizes the space-time evolution of the QGP. Considering 
that dimuonium is a heavy bound state, we use a classical 
Boltzmann-type transport equation to describe its evolu- 
tion. The distribution /(pt, y, xt, 77, T|b) as a function of 
transverse momentum pt and coordinate Xf , longitudinal 
rapidity y and space-time rapidity 77 and proper time t 
at fixed impact parameter b in a heavy ion collision is 
characterized by the equation 



cosh(?/ - 7?)-^ + - sinh(y -ri)S- 
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f = f3 



(4) 

with the dimuonium transverse velocity Vj = pt/Et and 
transverse energy Et — i/mj+pf. Because the life 
time of dimuonium is much longer than the life time 
of QGP, the decay of dimuonium is ignored. Moreover, 
since dimuonium does not participate in strong inter- 
actions, the generated dimuonium interacts with QGP 
only electromagnetically, the dissociation can then be ne- 
glected too. Therefore, we consider only the gain term 
/3(pt, y, xt, 77, r|b) on the right hand side of the transport 
equation, 

1 r (fk (fqi d^<i2 
2Et 
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where Eg = ym^-l-qf. Eg = ^ml+c^ and Eg = 

-(- are respectively quark, anti-quark and gluon 
energies, and /g, fq and fg are the thermal distribu- 
tions for quarks and gluons, fq — 1/ ^e'^i "^"/^ -f 1^ , 

fq ^ 1/ (e92"M/r + 1) and fg = 1/ (e'^''"^'/'^ - l). By 

using Bjorken's hydrodynamics [ll], the fluid velocity 
appeared in the distributions and the temperature T in 
the distributions and gluon mass are functions of coor- 
dinates (xt , rf) at fixed b and determined by the ideal 
hydrodynamic equations fl?'] 

drE + V-m = ~{E+p)It, 



3 



drMy + V • (AfyV) = -My/T " OyP, 

drR + V ■ {R^^) = -R/t 



(6) 



with the Lorentz factor 7 = — and definitions 

E = {e+p)j^—p, M — (e+p)7^v and R = jnaa functions 
of energy density e, pressure p and baryon density n of the 
medium. To close the hydrodynamic equations, we need 
a equation of state to describe the nature of the QGP. We 
take the resuh from the lattice QCD simulation with a 
phase transition temperature of deconfinement = 190 
MeV [1^. Since the colliding energy is so high in heavy 
ion collisions at RHIC and LHC, the net baryon density 
in the fireball is rather small [HI , we can simply set n = 
in the numerical calculations. The initial condition of 
the hydrodynamic equations is controlled by the nuclear 
geometry which determines the impact parameter and 
the colliding energy which governs the ratio of soft to 
hard contributions [11] . 

Since we did not consider the initial dimuonium pro- 
duction before the QGP formation and have neglected 
the loss term in the transport equation its analytic 
solution becomes simple, 

/ Pt,y,xt,77,T b = / dr — — 

Jr„ A(t') 

xe(r(Xt,i/,T'|b)-T,) (7) 

with the definitions 



Xt(r') = Xi-Vt[TCOsh(2/-r/)-r'A(r')], 
H{t') ~ y ~ arcsinh (r/T'sinh(y — r/)) , 

A(t') = ^l + (r/T')2sinh2(y-,7), 



(8) 



where the local temperature T(Xt, iJ, r'jb) as a func- 
tion of time and coordinates at fixed impact parameter 
is determined by the evolution of the medium dH), the 
step function indicates that the mechanism of gen- 
erating dimuonium discussed here can only take place in 
deconfined region, and the coordinate shifts Xj — > Xj and 
77 — >■ in the solution ([7|) reflect the leakage effect in the 
transverse and longitudinal directions. The time integra- 
tion is from the initial time tq to t. By integrating the 
distribution over the phase space, we obtain the dimuo- 
nium transverse momentum distribution at fixed impact 
parameter b. For the finally observed dimuonium distri- 
bution, r should be so chosen that it is not earlier than 
the end time r/ of the QGP. 
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where is an arbitrary time after Tf, since the momen- 
tum distribution is unchanged for all t > Tf. 



We can define the transverse energy distribution 
f{Et\h) = dN{h)/{2TTEtdEt). It is shown in FigH 
for central (b=0) Au+Au collisions at RHIC energy 
^Jsnn = 200 GeV and Pb-|-Pb collisions at LHC energy 
^Jsnn — 5.5 TeV. We have taken the effective coupling 
constant = 0.3 (corresponding to T /Tc — 1.5 — 2 [13]) 
and the dimuonium mass ~ 2m ^ = 211 MeV (ne- 
glecting the binding energy 1.4 KeV) in the numerical 
calculations. While the dimuoniums distribute in a wider 
region at LHC, they behave similarly at two energies. 
The result in the low Et region oi Et <\ GeV can be pa- 
rameterized as a thermal distribution f{Et) ^ e~^'^^'^'^' 
with a slope parameter Tg// = 195 MeV at RHIC and 240 
MeV at LHC. This indicates that the thermodynamic in- 
formation of the medium carried by quarks and gluons is 
partly inherited by the produced dimuoniums and can be 
used to signal the QGP formation in high energy nuclear 
collisions. 
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FIG. 2: The dimuonium transverse energy distribution 
dN / {2n EtdEt) in central Au-|-Au collisions at RHIC energy 
y^SNN = 200 GeV and Pb-|-Pb collisions at LHC energy 
y/s^ = 5.5 TeV. 




FIG. 3: The rescaled dimuonium number lO'^A^ as a function 
of participant number Np in Au-|-Au collisions at RHIC en- 
ergy ^ysNN ~ 200 GeV and Pb-(-Pb collisions at LHC energy 
ysTTjv = 5.5 TeV. 
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FigIS] shows the momentum integrated dimuonium 
yield N{h) = J f{pt\h)ptdpt as a function of the num- 
ber of participant nucleons Np in heavy ion cohisions at 
RHIC and LHC energies. The relation between Np and 
the impact parameter b can be easily determined by the 
nuclear geometry. With increasing centrality, the par- 
ticipant number increases, and the temperature, the life 
time and the space region of the formed QGP increase. 
As a result, the dimuonium yield goes up with centrality, 
due to the enhancement of the quark and gluon numbers. 
For central collisions with maximum Np, the dimuonium 
number is 1.3 x 10~^ at RHIC energy and becomes about 
3 times larger at LHC energy, as shown in Figl3l 

Let's now compare the two dimuonium production 
mechanisms in heavy ion collisions. One is through the 
pure electromagnetic channel, A1A2 — >■ AiA2{fJ,^ fi^), 
proposed by Ginzburg et al. !9|, and the other is inside 
the formed QGP, qq — >■ {^~^^~)g, discussed here. For 
impact parameter b > 2Ra, where Ra is the colliding 
nuclear radius, there is no QGP formed, the production 
is only through the electromagnetic channel. However, 
for h < 2Ra, the colliding nuclei are broken, and the 
assumption of the replacement of the perturbation pa- 
rameter a by Za with each photon exchange is no longer 



valid for the electromagnetic channel. For b <C 2Ra, the 
QGP is formed with a high temperature, long life time 
and large size, the process qq — > {fi'^^^)g becomes dom- 
inant. Since the pure electromagnetic channel is not re- 
lated to the fireball, the transverse energy distribution in 
this channel can not show the thermodynamic behavior. 

In summary, we investigated the dimuonium produc- 
tion in the QGP formed in relativistic heavy ion colli- 
sions. The dimuonium motion in the QGP is described 
by a transport equation with the gain term character- 
ized by the production process qq — >■ {^'^ ^~)g, and the 
space-time evolution of the plasma is controlled by ideal 
hydrodynamic equations. By solving the coupled trans- 
port and hydrodynamic equations for high energy nuclear 
collisions at RHIC and LHC energies, we found that while 
the electrodynamics dominated dimuonium yield is not 
high enough, the transverse energy distribution inherits 
the thermodynamic behavior of the hot medium and can 
be considered as an electromagnetic probe of the QGP. 
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